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ABSTRACT 

We report the results of a spectroscopic study of the Bootes III (BooIII) stellar 
overdensity carried out with the Hectospec multifiber spectrograph on the MMT tele- 
scope. Radial velocities have been measured for 193 BooIII candidate stars selected 
to have magnitudes and colors consistent with its upper main sequence and lower red 
giant branch, as well as a number of horizontal branch candidates. From 20 identified 
candidate BooIII members, we measure a systemic velocity of V & = 197.5 ±3.8 km 
s" 1 and a velocity dispersion of a Q = 14.0 ±3.2 km s" 1 . We use the somewhat large 
velocity dis persion and the i mplied highly radial orbit, along with morphological evi- 



dence from lGrillmairl (|2009[) and stellar abundances, to argue that BooIII is likely the 
first known object observed in a transitional state between being a bound dwarf galaxy 
and a completely unbound tidal stream. 
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namics - Galaxy: structure - Local Group 
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1. Introduction 



Over the last few years and thanks to the advent of large-area photometric sur veys such as the 



Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SPSS: lAbazajian et al. 



2009b . the substructured nature of the Milky Way halo has been clearly revealed. In r ecent years the 



number of known Galactic dwarf spheroidals (dSphs) ha s more than double d (e.g.. lWiUman et al. 



2005; Bel okurov et alJl200 6a. 2007; 



streams have been identified (e.g. 



Zucker et al 



Ibata et al. 



Yannv et al.ll2003l : Belokurov et alJl2006tt iGrillmaii 



2006; 



2001 



E 



Grillmair 



Maiewski et al. 



2003: 



2006a ), and numerous stellar 



Rocha-Pinto et al. 2004; 



2006 alibi 120091 : iGrillmair & Dionatosll2006n 



These discoveries lend significant support to currently favored cold dark matter cosmological mod- 
els, which predic t that spiral galaxies form hierarchically from t he merger and accretion of smaller 
subsystems (e.g.. bullock et aLlboOlUBullock & Johnstonll2005|) . 



One exciting aspect of these discoveries is that they have revealed a previously unknown 
population of "ultra-faint" dwarf galaxies which are claimed to contain the highest dark mat- 
ter fractions reported to date, with mass-to-light ratios abov e 1,000 (M/L) Q in some cases (e.g., 
Munoz et al.ll2006l; Isimon & Gehall2007l ; Istrigari et al.ll2008|) . but luminosities as low as L v ~ 300 



L (IMartin et al 



2007|) . A few of these satellites, Coma, Segue 1, Segue 2 and Bootes II, have 



properties closer to those of globular clusters than those of "classical" dSphs (e.g ., half-light radii 



betwe en 30 and 70 pc), and have been suggested to be satellites of satellites (|Belokurov et al 
20091) . 



In this Letter we present the results of a sp ectroscopic surv ey of another potential dwarf satel- 
lite, the Bootes III stellar overdensity (B00III; lGrillmairil2009L hereafter G09). This peculiar ob- 
j ect, located at a dis tance of 46 kpc and with V-band absolute magnitude of M v = -5.8 ±0.5 
(ICorrenti et al.ll2009l) . has a roughly power-law surface brightness profile and an irregular shape. 
Our study suggests that B00III is likely a dwarf satellite in a transitional state between a bound 
entity and an u nbound stellar stream. In this context, B00III resembles the objects described by 



Kroupal (|1997|) but likely with a high dark matter content. 



2. Spectroscopic Observations and Data Reduction 



Radial velocities (RVs) wer e measured from spe ctra obtained in 2009 February with the Hec- 
tospec multifiber spectrograph (|Fabricant et al.ll2005l) mounted at the f/5 focus on the 6.5 m MMT. 
Targets were selected from the SDSS stellar database to lie on or near the B00III color-magnitude 
sequence from IG09L in a magnitude range of 18. 5< g <22.5. We also i nclud ed candidate blue 
horizontal-branch (BHB) stars that are near the BHB locus in Figure 13 of lG09L bringing the total 
number of assigned targets in the ~1 degree Hectospec field of view to 227. Figure 1 shows the 
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g—i versus g color-magnitude diagram (CMD) for all stars within 0.5° of the BooIII center, as well 
as an old, metal-poor isochrone chosen to be similar to that of Ml 5 (the cluster which guided our 
target selection). 

We used the 600 gpm grating, centered at 5800A, to give a working wavelength range of 
4550-7050A at a dispersion of ~0.55 A pixel" 1 (2.1 A resolution). This region was selected to 
include the H/3, Mg triplet, Na D, and Ha spectral features. Queue-observed integrations totaled 3 
x 1800 s on UT 2009 February 3, and 6 x 1800 s on UT 2009 February 20, for a summed exposure 
of 4.5 hr on this field. Sky fibers (57 in all) were assigned throughout the field, distributed so that 
a number of them would fall within each of the two CCD chips of the Hectospec system. A total 
of 193 target spectra had star-like spectra (there were a number of background QSOs/AGN), with 
sufficient signal-to-noise ratio (S/N > 10) for RV measurement. 

The Hec tospec data were reduced using an external version of the SAO "SPECROAD" reduc- 
tion pipeline ( Mink et al. 2007 ) written by Juan Cabanela, called ESPECROAE0. The pipeline au- 



tomates many reduction steps, including bias-correction, flat- fielding, cosmic-ray rejection, fiber- 
to-fiber throughput adjustments, and sky subtraction. Wavelength calibration was performed man- 
ually using sets of three combined HeNeAr calibration lamp exposures from each night of observ- 
ing. 

We derived RVs using the IRAF task FXCOR to cross-correlate object spectra against two 
template G8V RV standards, each observed three times. We first correlated the standard spectra 
against each other, and found that to achieve agreement between our measurements and published 
RVs^l for these stars, it was necessary to shift all RVs measured against one of them by 6 km s" 1 . 
After doing so, the measured RVs for all six standard spectra agree with published values to within 
1.0 km s" 1 . To minimize the effects of noise for spectra with lower S/N, the cross-correlation was 
restricted to the regions around the Ha, Mg triplet, and H/3 lines. The standard deviation of six RV 
measurements for most objects (after applying the 6 km s" 1 shift discussed above) was less than 1 
km s" 1 , increasing to ~3 km s" 1 at g ~22. 

To estimate RV uncertainties, we separated our nine individual object exposures into three 
sets of three summed exposures, totaling 1.5 hr of exposure each. We then measured RVs in each 
of these three independent frames, and fit a curve to the run of standard deviation of the three 
measurements as a function of magnitude. We used this fit to the trend of <trv versus g for the 
brighter stars in our sample to assign an RV uncertainty for each star we observed based on its 
magnitude (extrapolating the curve for fainter stars having insufficient signal to measure an RV 



^http : //iparrizar .mnstate . edu/~ juan/ research/ESPECROAD/ index. php| 
3 From the Geneva Radial Velocity Standard Stars at |http : //obswww . unige . ch/~udry/std/| 
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from the 1.5 hr exposures). Velocity uncertainties vary from ovheiio ~ 3 km s 1 at g =18.5 to ~10 
km s" 1 at g = 21 (and higher at fainter magnitudes). 



3. Spectroscopic Results 

3.1. Bootes III Membership 

Figure 2 shows RVs for all stars in our sample with reliable measurements (S/N > 10). For 
comparison, we show (red dot-dashed curve) a scaled RV distribution from ~100 realizations of the 



Besancon Galactic model (|Robin et al.H20 03) along this line of sight, covering the same magnitude 
and color ranges as our targets. Since this field is at high Galactic latitude (b ~ 75.4°), blue stars 
at faint magnitudes are mostly Milky Way halo main sequence turnoff stars, with a large (~100 
km s" 1 ) velocity dispersion. This halo distribution is well matched by our data set. However, two 
significant excesses with respect to the expected halo distribution, at V he iio of -180 km s" 1 and ~200 
km s" 1 , are visible in the histogram. To determine whether one of these RV signatures is associated 
with BooIII we check the velocities of the six candidate BooIII BHB stars we have included from 



the "clump" in Figure 13 of lG09i (asterisk symbols in Figure 1). Interestingly, all six of these stars 
have velocities within the ~200 km s" 1 peak, with a mean value of Vh e u = 195.6 ± 7.2 km s" 1 and 
dispersion 13.7 km s" 1 . Since none of the BHB candidates have velocities within the negative RV 
peak, we identify this positive, high-velocity peak with the BooIII kinematical signature. 

Having thus identified BooIII's systemic velocity, we select other stars within 3-cr of this BHB 
velocity peak as candidate BooIII members. These stars are highlighted in the g — i versus g CMD 
in Figure 1. The good agreement between the isochrone for an old, metal-poor ([Fe/H] = -2.30) 
population and the BooIII RV candidates suggests that we are indeed seeing the red giant branch 
of an old, metal-poor, comoving population of stars. 

Cross-correlation of BHB spectra in our sample against the two available G-dwarf standard 
spectra results in RVs with higher uncertainties than we derive for the later-type stars in our sample. 
This is due to the mismatch in spectral type and spectral line profiles between the BHB candidates 
and the RV standard stars. Because of this, the scatter in our BHB RVs is likely artificially high, 
potentially leading to higher field-star contamination in our 3-cr velocity selection. In order to reach 
a cleaner sample of BooIII members, we further exclude stars that are more than 0.25 magnitudes 
from the isochrone shown in Figure 1 . This leaves our final sample of 20 BooIII candidate members 
(14 if the BHB stars are excluded), which is represented in Figures 1, 3, and 4 with large filled 
symbols and asterisks. 
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Fig. 1. — SDSS color-magnitude diagram of a 30 ar cmin radiu s abou t the nominal BooIII cen- 
ter, drawn from the same data points as Fig ure 13 of iGrillmain (|2009l) . The solid line is an old, 
metal-poor ([Fe/H] = -2.30) isochrone from lGirardi et all J2000). shifted to a distance of 46.0 kpc. 
Large symbols represent all stars within the 3-cr RV selection centered at ~197.5 km s" 1 ; asterisks 
represent the six BHB candidates that were in our spectroscopic sample, triangles are RV members 
whose spectra had S/N > 20, and squares those having S/N < 20. Filled symbols and asterisks 
are candidate RV members that are also within 0.25 magnitudes of the metal-poor isochrone in the 
CMD; these constitute our final BooIII sample. 
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Fig. 2. — Heliocentric radial velocities (solid histogram) for all stars in the MMT+Hectospec con- 
figuration with reliable measurements. For comparison we show scaled output from the Besancon 
model as the red dot-dashed histogram. The model data consist of stars within 0.2 < g — r < 0.6, 
and 18.6 < g < 22.0, which are predominantly Milky Way halo turnoff stars in the same region 
of the CMD as our targets. Relative to the model, there are two peaks of "excess" counts: one at 
roughly - 1 80 km s" 1 , and a second at V~200 km s" 1 . We identify this latter peak as associated with 
BooIII. The blue dashed curve is a Gaussian with mean velocity V Q = 197.5 km s" 1 and a width of 
cr = 14.0 km s" 1 , matching our final derived values for BooIII. 
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3.2. Metallicity 



Stellar atmospheric parameters for all observed stars other t han the BHB candid ates are de- 
termined using low-resolution Lick spectral indices (as defined bv lWorthev et al.ll 19941). Measured 
indices for all input spectra are placed onto the standard Lick index system of ISchiavonl (|2007|) 
using offsets derived from standard spectra observed with the same telescope/instrument setup. 
We derive a suitable estimate of the effective temperature, r eff , of each star based on its (g - r) 
photometric color. Relatively accurate dwarf/giant separation is then performed using the gravity- 
sensitive Mg b feature at ~ 5170 A , with stars sorting well by gravity in the Mg b-H/3 index plane 
(H/9 is primarily s ensitive to r eff ). The metallicity, [Fe/H], for each star is determined using two 
methods; from the ISchiavonl (|2007r) spectral library we have derived separate polynomial surface 
fits of (1) [Fe/H] to the Lick i ron indices and stellar color, and (2) [Fe/H] as a function of the iron 
indices and H/j (sim ilar to the lFriellll987l interpolation between lines of constant metallicity, or the 
Gorgas et al.lll993l fitting functions). Thus, moderately accurate estimates of [Fe/H] derive from 
the combination of all six Lick iron indices falling within the observed spectral bandpass as well 
as the stellar color and H/9 index. 

Measured metallicities as a function of Vheiio are shown in Figure 3. We estimate typical 
uncertaint i es in [Fe/H] of ~0.3-0.5 dex, depending on the magnitude of the star. According to 
Grillmairj (|2009), the matche d-filter signa l of BooIII is strongest for the most metal-poor filter 
applied (M15: [Fe/Hl =-2.26: iHarrisI 1 1 9961) . Furthermore, the M15 ridgeline provides the best fit 
to the CMD locus of BooIII ( G09i Figure 12). We thus expect the majority of BooIII members to 
be metal-poor ([Fe/H] < -2.0). Examining the stars within the 3-a RV selection in Figure 3, we 
note that a number of candidates are indeed rather metal-poor. Furthermore, both of the stars with 
S/N > 20 spectra that were excluded from the final sample as likely foreground contaminants (i.e. 
open triangles in Figures 1,3, and 4) have measured [Fe/H] > -1.8. 

From the individual stellar metallicities of our final sample (excluding the BHB stars) we 
derive a maximum likelihood estimation of [Fe/H] « -2.1 ±0.2 for BooIII. Since these individual 
values are derived from low S/N spectra, we shift the BooIII candidate spectra to rest wavelength 
and stack them to produce a higher S/N spectrum. Doing so, we derive [Fe/H] « -2.0, which is 
consistent with all other me asurements for BooI II. At the integrated V-band absolute magnitude 



of BooIII, My = -5.8 ± 0.5 (Correnti et a 



2009), our measured metallicity is typical for a dwarf 



galaxy (see Figure 11 of Simon & Gehal 2007|) . and more metal-poor than all Galactic globular 



clusters at similar luminosity. As is obvious from Figure 3, there is a fairly large metallicity spread 
in the BooIII candidate stars (o^e/H] ~ 0-6). This may be due to foreground contamination of our 
sample, though as seen in Figure 2, very few Galactic stars are expected at these velocities. This 
spread in [Fe/H] is comparable to (though slightly higher than) those found by ISimon & Geha 
(|2007|) for a number of other ultrafaint dwarfs (see their Table 4), and may reflect the presence of 
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multiple generations of star formation among the BooIII stellar population. 



3.3. Velocity Dispersion and Mass Estimate 



With our final sample of BooIII candidates, w e use the maximum likelihood method (|Pryor & Meylan 



19931 : lHargreaves et al.ll 19941 : iKleyna et al.ll2002h to derive the systemic velocity and velocity dis- 
persion of Bootes III. For the full sample (20 stars, including the 6 BHB stars), we find V = 
197.5 ±3.8 km s" 1 and a velocity dispersion of a D = 14.0 ±3.2 km s" 1 . Removing the BHB stars 
(leaving 14 total), the result is virtually unchanged: V = 198.2 ±4.5 km s" 1 , and a Q = 14. 1 dz 3.7 
km s" 1 - we thus consider only the full sample henceforth. At the position and distance of BooIII, 
this Vq converts to an unusually high Galactic standard of rest velocity of V G sr = 238.8 ± 3.8 km 
s -1 which implies that BooIII is on a very radial (and thus potentially destructive ) orbi t. This idea 



of a high-latitude object on a rather eccentric orbit agrees with the suggestion by lG09l that BooIII 
may be extended along the line of sight (i.e., along its orbit). 

Our measured dispersion of 14 km s" 1 for BooIII is the highest measured line-of-sight velocity 
dispersion of any dSph, even higher than the central dispersion of Sagittar ius, the most massive 



of the Galactic dSph galaxies. It has been shown in previo us studies (e.g.. IPiatek & Pryorlll995 



Oh. Lin. & Aarsethll 19951 : iKroupall 19971 : iMunoz et al 



2008) that even if a satellite is experiencing 



tidal disruption, its cr remains unaffected by tides unless the system is at or near the point of 
complete destruction (as we believe is the case for BooIII). In this context, it is possible that even 
in the case o f Sagittarius t he cen tral dispersion reflects the equilibrium value. The same studies, 
in particular Munoz et al.l (|2008l) . show that when a satellite is on the brink of being completely 
destroyed by the action of tides, its central velocity dispersion inflates in some cases by more than 
a factor of 2. If BooIII was originally a system similar to the classical dSphs, with a half-light 
radius of ~ 300 pc and a velocity dispersion of ~ 6-7 km s" 1 , then its current measured value of 
14 km s" 1 is consistent with a scenario wherein BooIII is near complete destruction. 

In order to e stimate the dynamical mass of BooIII we use a method recently developed by 
Wolf et all i2009h which relies solely on the central velocity dispersion and the half-light radius 
of the system, and assumes neither that the density distribution of the satellite follows a King-like 
profile nor that mass follows light. We note that this does assume we are observing a bound system 
in dynamical equilibrium, which we argue is not the case for Bootes III. For illustrative purposes, 
however, we carry out the exercise. The expression derived by Wolf et al. depends only on direct 

3o LOS r l/2 4o-£ os R eff 



observables, and is independent of the velocity anisotropy: M\j% 



Here, <7 LO s 



G — G 

is the line-of-sight velocity dispersion, r\n is the three-dimensional deprojected half-light radius, 
and R e ff is the two-dimensional projected half-light radius (using their derived approximation R e f[ ~ 
(3/4)n/ 2 ). Substituting for the constants, this can be expressed as M(< n/ 2 )~9.3 x 10 4 |§f^(f^y|) 2 
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Fig. 3. — Le/j panel: heliocentric radial velocities as a function of metallicity ([Fe/H]) for all stars 
in our observed sample having spectra with S/N > 20. Right panel: RVs vs. SDSS g magnitude 
for all stars with spectra. Symbols for both panels are as in Figure 1. As in previous plots, open 
symbols are those stars that were excluded from the final sample of B00III candidates. The hori- 
zontal lines in both panels represent the mean (solid line) and 3 a (dot-dashed line) RV selection 
criteria based on the BHB stars. 
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M Q . Taking our measured velocity dispersion a = 14.0km s we find M(< ri/ 2 ) c± 1.82 x lO 7 -^^ 
M Q , where we've left R e ff explicit to allow for the as-yet unmeasured surface brightness profile of 
BooIII. 



This derived mass is rather typical for a dSph system - IStrigari et al.l (120081) showed that 
Local Group dwarf galaxies all share a common mass scale of the order of ~ 10 7 M^ measured 



within the inner 300 pc of these satellites (as had earlier been suggested by Mated 1 19981) . Note 
that the above mass estimate for BooIII would be three times higher if we take i? e ff =300 pc, 
although the true R e ff is likely to be smaller (|G09I) . Adopting M v = -5.8 gives a total luminosity of 
L = 1.8 x 10 4 L Q . Using this value along with our mass estimate, we derive a total mass-to-light 
ratio of (M/L) v ~ 1000-j^ (M/L) V)Q . A mass-to-light ratio of ^ 1000 is similar to those measured 



for the ultrafaint dSphs of similar M v (e.g. iSimon & Gehal l2007. Figure 15). 



4. Discussion 



Despite the seemingly "normal" mass derived for BooIII we argue that a number of pieces of 
evidence support the idea that BooIII is not in dynamical equilibrium and it is likely in a transitional 
state betwee n a bound syst em and a completely unbound stellar stream - similar to the objects 
described bv lKroupaldl997l) but with a high dark matter content. As we have already suggested, the 
fact that BooIII has the highest m easured line-of-sigh t velocity dispersion of any dSph is consistent 
with the observation (from, e.g- JMunoz et al.ll2008|) that a satellite in the throes of complete tidal 
disruption can have its velocity dispersion inflated by more than a factor of 2. Tidal disruption 
seems quite plausible, since the large measured Galactocentric velocity for BooIII places it in a 
highly elongated, and thus potentially destructive, orbit. 

Morphological evidence further supports the suggestion that Bootes III is a d issolv ing dwarf 
galaxy in the final stages of ti dal destruction. The density distribution derived by lG09l (and con- 
firmed by ICorrenti et al.ll2009|) is quite unusual. Most dSphs, with perhaps the exception of a few 
of the new ultrafaint ones which are suspected to be undergoing tidal stripping, have density pro- 
files that are well d escribed by King-lik e or Plummer-like profiles. In contrast, Figure 1 1 of 1g09 



(and discussion by ICorrenti et al. 2009 ) shows that BooIII is better described by a power law of 
index -1. This observational result is also consistent with the idea of a satellite on the brink of 



destruction. Several studies (e.g.. IJohnston et alJI 1999c iMunoz et al.ll2008c IPenarrubia et al.ll2008l) 
have shown that a satellite which started with a King-like density distribution develops a power 
law-like distribution in the outer parts of its luminous component as it undergoes tidal stripping 
(due to escaping stars at large radii). Eventually, and as the satellite becomes nearly unbound, the 
King-like core also turns into a power-law like distribution. Furthermore, Grillmair (2009; Fig. 
10) showed that BooIII has irregular stellar surface density contours and is extended along the 



210.0 



209.5 209.0 
RA (degrees, J2000) 



Fig. 4. — Spatial distribution of all observed spectroscopic targets. Symbols as in Figure 3. Con- 
tours represent filtered star counts, and the large plus sign denotes the BooIII center, both from the 



G09| study. 
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east-west direction. In Figur e 4, we show the spatial distribution of our spectroscopic targets over- 



laid on contours of the |G09J filtered star counts. Though our single MMT+Hectospec setup covers 
only a fraction of the rather large BooIII stellar overdensity, the identified RV members seem to 
follow the BooIII contours. The Styx stellar stream also passes through the same region of the sky 
where BooIII is located, and at nearly the same distance (IG09 ). BooIII may be associated with this 
stream, and indeed may even be its progenitor. However, further kinematical information along 
the Styx stream will be necessary to confirm the association. 

We thank the anonymous referee for helpful suggestions, and gratefully acknowledge support 
by NSF grant AST-0807945 and NASA/JPL contract 1228235. 
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